ftQ  (eft  ft# 


This  document  has  been  approved 
for  public  release  and  sale;  Its 
distribution  Is  unlimited. 


AD 


TECHNICAL  REPORT 
69-51-OSD 


APPROXIMATE  ANALYSIS  OF  A  FLAT,  CIRCULAR  PARACHUTE 
IN  STEADY  DESCENT 


by 


Edward  W.  Ross,  Jr. 


December  19 68 


Office  of  the  Scientific  Director 
U.S.  ARMY  NATICK  LABORATORIES 
Natick,  Massachusetts 


Best  Available  Copy 


I 


FOREWORD 


Despite  the  fact  that  parachutes  have  been  in  use  for  many  years, 
no  completely  satisfactory  theory  exists  for  the  prediction  of  stresses 
in  a  canopy.  This  hinders  good  design.  This  report  describes  a 
mathematical  model  for  predicting  canopy  shape  and  stresses  in  steady 
descent  which  is  superior  to  those  now  in  use.  This  model  is  embodied 
in  a  computer  program  that  can  be  used  for  design  studies  if  desired. 
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ABSTRACT 


A  theory  is  presented  for  the  stress  analysis  of  a  flat, 
circular  parachute  in  steady,  vertical  descent.  Unlike 
previous  treatments  of  the  problem,  this  theory  does  not 
assume  that  the  shape  is  known.  Instead5the  theory  presents 
relations  between  the  pressure  distribution  in  the  opened 
condition  and  the  shape,  drag  and  stresses  in  lines  and 
fabric.  The  theory  results  in  a  non-linear  third  order 
system  of  ordinary  differential  equations  with  boundary 
conditions  at  both  vent  and  skirt.  This  system  was  solved 
by  a  computer  program  based  on  the  Runge-Kutta  method  of 
numerical  integration.  The  results  are  in  fairly  good 
agreement  with  measurements  on  parachutes.  The  computer 
program  can  be  used  for  studies  of  effects  of  design  changes 
on  shape,  drag  and  stress,  and  the  results  of  a  small  study 
of  this  sort  are  included. 


1.  INTRODUCTION 


In  recent  years  repeated  attempts  have  been  made  to 
analyze  the  behavior  of  parachutes  and  devise  formulas 
suitable  for  their  design.  Much  progress  has  been  made, 
but  the  state  of  knowledge  is  still  not  wholly  satisfactory. 
The  best  analyses  currently  available,  Heinrich  and  Jamison* 
snd  Topping,  Marketos  and  Costakos  rely  e:  knowing  the 
deformed  shape  as  well  as  the  pressure  distribution  before 
the  stresses  are  analyzed  even  in  steady  descent.  Present 
knowledge  about  large  deflections  of  structures  suggests 
that  it  should  be  possible  to  calculate  the  stresses  and 
deformed  sh-os  concur  rent  ly,  although  some  assumptions  must 
still  be  made  about  the  pressure  distribution  in  the 
deformed  state.  We  shall  undertake  to  do  this  in  the 
present  paper. 

To  be  specific,  we  analyze  a  flat  circular  parachute  in 
steady,  vertical  descent  by  an  approximate  theory  for  large 
elastic  deflections.  This  theory  is  like  (but  goes  beyond) 
that  of  Heinrich  and  Jamison*.  The  parachute  is  regarded  as 
a  completely  flexible  structure,  i.  e  none  of  its  elements 
(cloth,  cords,  reinforcing  tape,  etc.)  have  any  bending 
stiffness.  The  resulting  analysis  resembles  Che  large- 
deflection  (non-linear)  membrane  version  of  thin  shell  theory 
but  differs  from  it  because  of  the  important  pare  played  by 
the  cords.  A  number  of  assumptions  are  made,  of  which  the 
most  important  are  listed  below: 
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(1)  The  strains  In  the  fabric  and  cords  are  small  even 
though  the  deflections  are  large. 

(11)  The  fabric  is  supposed  to  possess  no  resistance 
to  stresses  in  the  meridional  (radial)  direction,  all  forces 
in  this  direction  being  resisted  solely  by  the  cords. 

(ill)  Corresponding  points  on  each  gore  experience 
identical  stresses  and  deformations,  and  the  same  is  true 
for  each  cord,  i.  e the  deformation  is,  loosely  speaking, 
axi ally-axisymraetri c. 

(iv)  The  circumferential  radius  of  curvature  of  the  gore3 
is  everywhere  much  smaller  than  the  meridional  radius  of 
curvature . 

(v)  For  each  gore  points  on  a  circular  arc  about  the  axis 
of  the  undeformed  parachute  lie  in  a  plane  perpendicular  to 
the  cords  after  deformation. 

The  analysis,  which  is  described  in  the  next  Section, 
leads  to  a  non-linear  third-order  system  of  ordinary  differential 
equations,  which  cannot  be  integrated  in  closed  form.  Therefore 
a  computer  program  was  written  for  the  solution  of  the  system, 
based  on  the  Runge-Kutta  method  of  numerical  integration.  The 
program  Is  automatic  enough  so  that  it  Is  convenient  to  study 
the  effect  on  the  deformation* ,  stresses,  drag  and  parachute 
weight  of  any  two,  arbitrarily  chosen,  design  parameters,  such 
as  the  vertical  velocity,  number  of  gores  or  elastic 
moduli  of  various  structural  elements.  A  description 
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of  the  program  is  given  in  Section  3  ,  and  examples  of  the 

results  obtained  are  displayed  in  Section  4.  Section  5 
contains  a  discussion  of  the  analysis  and  results. 
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2.  ANALYSIS 


In  analyzing  the  parachute  we  first  treat  a  generic  one 
of  the  gores,  which  is  assumed  to  be  a  fiat  sector  of  a  circle 
in  the  undeformed  state.  If  G  is  the  number  of  gores,  then  the 
sectorial  angle  in  both  the  undeformed  and  deformed  states  is 

o  -  itr/G  (1) 

We  examine  first  the  kinematics  ,  then  the  statics  and  finally 

the  material  behavior  of  the  fabric. 

Kinematics:  Let  A  and  B  denote  the  material  points 

where  a  generic  circular  arc  in  the  undeformed  state  meets 

the  cords  that  border  the  gore.  In  the  undeformed  state  this 

arc,  AB ,  has  length,  R  a;  see  Figure  1.  We  assume  that  after 

deformation  the  arc  AB  forma  a  plane  circular  arc,  the  plane 

of  the  arc  being  perpendicular  to  the  deformed  cords.  Figure  2. 

The  points  A  and  B  are  now  at  a  distance  r  from  the  axis  and 

are  separated  by  a  distance  ra.  If  r  is  the  radius  of  this 

AB 

circular  arc  and  2g  is  the  sectorial  angle,  then  ws:  see  from 

Figure  3  that,  provided  0^b<u/2, 

r  sing  -  (l/2)ra  (2) 

AB 

2gr  »  Deformed  length  of  arc  AB  (3) 

AB 

!f  b  •  v/2,  we  may  have  a  situation  where  adjacent  gores  are 
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Figure  1: 


Undefonned  parachute  go 
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pressed  against  each  other  with  a  contact  length  a*  as  shown 
in  Figure  4. 

In  this  case 

rAB  *  C:/2)r«  (4) 

Deformed  length  of  AB  ■  2o'  +  tr 

-  2cj  '  +  ( tt/2  )  ra  (5) 

The  undeformed  length  of  AB  is  Ro,  so  we  may  write 
formulas  for  the  circumferential  fabric  strain,  which  is 

assumed  small, 

« 

Yf  ■  change  in  length/original  length 

Yf  -  (2BrAB/Ro)  -  1  if  B  <  it  /  2  (6) 

*  {[2o,  +  (tr/2)ra]/Ra>  -  \  if  B  •  tt/2  (7) 

Statics:  If  Nf  is  the  circumferential  tension  in  the 

gore  fabric,  with  dimensions  (lbs/ft)  and  p  is  the  pressure 
difference  between  the  Inside  and  outside  of  the  surface, 
then  equilibrium  requires 

Sf  '  "rAB  (8) 

The  forces  exerted  by  the  fabric  on  the  cords, 
per  unit  length  of  cord  are  shown  in  Figure  5.  The 
forces  in  the  circumferential  and  normal  directions, 

Nq  and  Nn f  ire  found  from  conditions  of  equilibrium, 

KS  “  Nf  coa  6  (9) 


Nn  -  2Nf  8in6 


(10) 
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Figure  5: 


Forces  transmitted  from  gore  fabric  to  cords 
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If  we  use  (2)  and  (8),  we  obtain  for  6  i1  0 


N0  “  ( 1 / 2 ) p r a  cot  8 

(11) 

N 

”n  “  ora 

(12) 

These  formulas  are  valid  both  for  B<*/2  and  8  *  ir/2,  and 
in  the  latter  case  they  imply 


Material  Behavior:  We  assume  that  the  fabric  obeys  a 
linear  elastic  law, 

Nf  "  EfYf  (13) 

where  is  the  elastic  constant  of  the  fabric  when  stressed 
in  the  circumferential  directicn.  If  8<*/2,  we  may  use  (2),  (6) 
and  (s)  to  rewrite  this  relation  in  the  form 

r/R  -  B"l{sinS+(l/2) (p/Ef)ra}  (14) 

This  can  also  be  written,  with  the  aid  of  (11)  S3 
r/R  -  B-’sin  8 { l+(No/Ec ) se c  B), 
or 

N  /E  -  cos  6(8  esc  8  (r/R)  -1*  (15) 

0  I 

We  see  from  (15)  that  in  the  limit  aa  fi  •+  0  we  obtain 
Nt//Ef  -*■  (r/R)  -1. 

This  can  be  interpreted  to  mean  that,  when  the  fabric  is  stretched 
flat  between  the  cords,  i.  e.,  when  there  is  negligible  bulging 
of  the  gores,  we  recover  the  linear  elastic  relation  between 
circumferential  stress  resultant  and  circumferential  atrain. 
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In  the  calculations  we  shall  adopt  (15)  together  with 
formulas  equivalent  to  (14)  and  (13) 

sin  8  -  8(r/F^f(l/2)  (p/Ef  )ra  -  0  (16) 

N  -  E  (8csc8(r/R)  -1}  (17) 

as  expressing  the  material  behavior  of  the  fabric  when 

B<*/2. 

When  8  ■  ir/2,  the  elastic  law  (13)  implies 
( 1/2) (p/E^) ra  *  { [  2o  '  +  (  ir/2)  ra]  /Ra  }  -1 
or 

o’  -  (1/2)Rq(1  -  ( ir/2)  (r/R)  +  (l/2)(p/E  )ro)  (18) 

together  with 

N  -  0  (19) 

0 

Nf  -  (l/2)pra  (20) 

We  have  now  completed  the  portion  of  the  analysis  that 
deafevlth  the  fabric. 

In  analyzing  the  fabric  we  have  neglected  the  meridional 
stresses  on  the  assumption  that  stresses  of  this  type  are 
borne  entirely  by  the  cords.  We  must  now,  therefore,  analyze 
the  statics  and  material  behavior  of  the  cords. 

Statics  of  Cords:  If  N£  is  the  tension  force  in  the 
cords,  measured  in  lbs,  then  we  can  write  down  the  two 
equations  of  force  equilibrium  in  the  direction  tangential 
and  normal  to  the  cords  with  the  aid  of  Figures  6  and  7. 

These  are 

dN  / ds  *  2N  sin(a/2)cos$  (21) 

c  e 

N  d$/ds  “  N  -2N  sin(a/2)sin4> 
c  n  u 

-  pra-2N0sin(a/2)sin$  (22) 
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Figure  6:  Sketch  showing  that  the  resultant  of  the  N 
forces  oh  a  cord  is  a  radial  irrwazxl  force.® 
2>;a  sin(a/2)~per"uhit  of  cord  length_ 
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Figure  7:  Forces  acting  on  an  element  of  cord  with  length  ds 
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A  partial  check  on  the  validity  of  these  equations  is 

provided  by  studying  their  behavior  as  or*  0 .  If  we  set  N  *  ra  N  , 

c  s 

these  should  reduce  to  the  equations  of  equilibrium  of  membrane 
shell  theory, 

d( rN  ) / dr  ■  N  (23) 

s  e 

d(sin$)/dr  +  Ng(sin  o/r)  -  d.  (24) 

If  we  use  the  fact  that  (see  Figure  8) 

dr  *  ds  cos  $,  (25) 

and  observe  that 

(2  /  a)sin(a/2)-+l 

as  a-*-0  ,  we  see  that  equations  (21)  and  (22)  reduce  to  the 
proper  limit  as  a-*0 . 

Material  Behavior  of  Cords:  The  strain  in  the  cords, 

Y3  ,  is  assumed  to  be  small  and  we  infer  from  Figure  9  that 

Y  ■  (ds-dR)/dR  z  (dr/ dR)  se  c<*>-  1 
s 

The  cords  are  assumed  to  obey  a  linear  elastic  relation 

N  =  E  y  -  E  { (dr/dRJsec'J’-l  >.  (26) 

c  c  s  c 

This  completes  the  analysis  of  the  cords. 

To  put  these  equations  into  systematic  form  we 
observe  that 

dr/dR  *  { 1+ (N C/E  )  }  cos  $  ■  f  cos  $  (27) 

whe  re 

f  -  l+(Nc/Ec). 

Since 

d/ds- (dr/ds ) (d/dr)»(dr/ds) (dR/dr) (d/dR) 
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we  find  from  (25)  and  (27)  that 
d/da  -  (1/f)  (d/dR). 

Thua  the  equations  of  equilibrium  and  elasticity  for  the 
corda  form  the  ayatem 

dr/dR  “  f  coa $ 

d*/dR  -  (f/Nfi) (pr a-2N 0a in ( a/2 ) a  in ^ } 
dNc/dR  -  2fN0ain(a/2)cos* 

where 

f  -  l+(Hc/Ec) 

and  N0  is  determined  as  followa.  If 
l+(r/2){(po/Ef)-(w/R) }  <0 

then 


cos6  ( (r/R)  8c«  cg-1 } 

where  B  is  found  by  solving  the  transcendental  equation 
ainB  -  6(r/R)  +  (pro/2Ef ).  -  0 
If 


l+fr/2){(po/Ef)  -  ( ir/R)  } >0 

then 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


N6  *  0  (34) 

a*  -  (Ra/2) (l+(pra/2Ef)  -  (*r/2R) } 

This  la  a  third  order,  non-linear  differential 
equation  ayatem  in  the  three  unknowns  r,  $  and  Nc,  and 
therefore  we  shall  need  three  edge  conditions  in  order 
to  have  a  unique  solution.  At  the  vent,  r  -  r^  or  R  ■  Rj, 
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the  cords  continue  across  to  the  opposite  side  of  the  vent 
through  the  axis  of  symmetry  (see  Figure  10).  The  elastic 
law  for  the  cords  implies  that 


r  -  R.+  {(N  ( R  )  /E  }  (35) 

i  1  c  i  c 

We  shall  take  for  another  condition  at  the  vent 

$  -  0.  (36) 

These  conditions  are  satisfied  with  good  accuracy  in  a  flat 
circular  canopy  although  in  other  kinds  of  parachutes, 
notably  those  with  the  pull-down  vent,  more  complicated 
conditions  must  be  imposed. 

The  third  condition  is  at  the  skirt,  r  -  r  or  R  ■  R  . 

o  o 

From  the  geometry,  see  Figure  11,  we  obtain 

r  -  L  cos  (rr-$_)  ■  r  +L  cos  4  *  0  (37) 

o  °  o  o 

where  L  is  the  deformed  length  of  the  suspension  lines 
between  the  skirt  and  the  load  (or  cor  fluence  point). 

If  Lq  is  the  undeformed  length  of  the  suspension  lin^s  , 
and  Y^is  the  strain  in  the  suspension  lines  (which  we  do 
not  assume  small),  then 

1/2 

L  -  L  ( 1+2 Y. )  ' 
o  u 
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The  tension  fcrce  In  the  lines  Is  simply  N  (Rq)  ,  the  tension 

force  in  the  cords  at  the  skirt,  hence,  if  E  is  the  elastic 

L 

modulus  (lbs)  of  the  lines. 


rt  ■  VV/Et 

<■  ■  i0d+<2vv/eJ1/2 

Then  from  (37)  we  find  the  condition  which  has  to  be  satisfied 
at  the  skirt. 


r (R  )+L  {  1+f 2N  (R  ) /Er ] } 1 1 2  cos* (R  )  -  0  (38) 

O  q  C  O  L  o 

We  have  therefore  reduced  the  problem  of  finding  the 

deformed  shape  and  stresses  in  the  parachute  to  that  of 

solving  the  system  of  three,  non-linear,  ordinary  differential 

equations  (28)  -  (30)  for  r(R)  and  ♦(R)  and  The 

quantities  f  and  ^  which  occur  in  the  right  sides  of  (28)  - 

(30) ,  are  determined  as  functions  of  Nc,  r  and  R  by  means  of 

(31)  •*  (34).  The  edge  conditions  on  the  differential  equation 


system  are  (35),  (36)  and  (38). 

When  this  system  has  been  solved,  the  drag,  or  weight  of 
the  load,  can  be  found  from 


D  -  GN  (R  ) s in$ (R  )  (39) 

c  °  o 

Also  the  deformed  shape  is  found  by  calculating  the  cord 

profile,  r(R)  and  Z(R),  and  the  gore  centerline  profile, 

r  (R)  and  Z  (R)  ,  see  Figure  12.  In  this  calculation  r(R) 

8  g 

is  of  course  known  directly  from  the  differential  equation 
solution,  but  Z(R)  must  be  found  by  integrating 
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gore  profile 


Figure  12:  Profiles  of  deformed  cord  and  gore  centerline 
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d7.  •  dr  tan  4  *  (dr/dR)  dR  tan  $ 


and  so  vc  find  with  the  aid  of  (28) 


Z  -If  (R’ )sin4(R' )dR* 


'R '  -  0 


Z  «\  { 1+(N  (R ' ) /E  JslntCRMdR' 

1  c  c 


R’  -  0 


The  gore  centerline  profile  la  then  found,  as  In  Figure  12, 
r  -  r+6aln$  (41) 

O 

Z  -  Z-6cos$  (42) 

g 

where  6,  the  depth  of  centerline  bulging.  Is  given  by 

different  formulas  for  a'  -  0  and  o’> 0.  If  o’  -  0, 

8<ir/2  and  (tee  Figure  3) 

6  -  r  (l-co«8)  ■  (1/2)  ra(l-cos8)/sin8 
AB 

-  (1/2) rat*n(S/2)  (43) 

If  a’>0,  8  •  if/2,  we  find  from  Figure  4 

6  -  (1/2)  ra+a'.  (44) 

We  shall  now  re-phrase  this  entire  system  of  equations  in 
dlaenslonless  fora  by  Introducing  the  new  definitions 


t  -  »Ri  , 

■e  ’  W 


r  -  yR£ 

N9  *  Efnc’ 


L  -  JR  ,  R  •  x  R  . 

0  1  o  o  i 

N  -  E  n  ,  7  •  ?  o 

f  ft'  L  ~  .  » 


P  *  «*2Ef/Rt’  °f  "  EfRi/Ec*  ut  -  E  /E 

“  L  c 


D  -  E  w, 

C 


oRj  , 


8  g  1  (45) 

’  rs’y*Ri’ 

6  -  AR. 


The  equations  become 


dy/dx  *  f  cos  $ 

(46) 

d$/dx  “  (2U,f/n  ){qay-n„  s  in  (  a/2  )  s  i  n  <$> } 
f  s  c 

(47) 

dn_/dx  ■  2fU,n  s in ( a/2 ) cos t 
s  f  c 

(48) 

f  -  1+n 

s 

(49) 

If 

l+a ay- [ ny / (2  x)  1  <0 

then 

n  ■  cos  8  {  8  (y /x)  cs  c8  -1} 
c 

(50) 

where  6 

satisfies 

sin  8  -  8y/x-qay  «  0 

(51) 

and  if 

1+q ay-  (iry/2x)  >0 

then 

n£  -  0,  8"x/2 

(52) 

a"  ( 1/2) xa { 1+q ay- ( ny /2x)  } 

(53) 

the  edge 

conditions  are 

ns(i)  -  y ( l) - 1 

(54) 

♦  (1)  -  0 

(55) 

1  /  2 

y  (xD)+J  {i+[2ng  (xq) /UL  ]  }  cos$(xQ)  -  0 

(56) 

The  dimensionless  load  or  drag  is 

W  *Gn  (x  )sin4>(x  )  , 

S  *■'  0 

(57) 
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Z(x) 


and 


“J"  { 1+n  (x’)}sin#(x')dx' 

(58) 

y+As  in  <J> 

(59) 

z-Acos  $ 

(60) 

(  1 /2  )y  at  an  (  6 /2  )  if  6<ir/2 

(61) 

o+(l/2)ya  if  B  »  n/2 

(62) 

In  discussing  the  pressure  and  stress  distributions  it  is 
also  convenient  to  define  dimensionless  stress,  drag  and  pressure 
coefficients  in  a  manner  different  from  the  dimensionless 
quantities,  previously  defined*.  We  shall  take 

2 

C  ■  p  /  {(1/2)  pv  }  ■  pressure  coefficient 
P 

C  ■  D/ { ( 1 /2) pV 2nR  2}  -  drag  coefficient 
Do  o 

n  2 

C  «  N  /{  (l/2)pV~  irR  }  ■  cord  stress  coefficient 

c  c  o 

2  2 

C,  -  NfR  /{(l/2)pV  irR  }  -  fabric  stress  coefficient 

t  r  o  o 


If  the  foregoing  analysis  is  correct,  we  must  satisfy  the 
3rd  order  differential  equation  system  (46)  -  (49)  together  with 
the  auxiliary  equation  (50)  -  (53)  and  the  edge  condition  (54)  - 
(56).  This  system  was  taken  as  the  basis  for  the  computer  program. 

The  numerical  analysis  underlying  the  computer  program  is  described  in  tbe 
next  Section. 

♦Notice  that  because  of  (13),  (26)  and  (35)  the  dimensionless  stresses 
and  n  ,  defined  in  (45),  are  also  respectively  just  the  strains 

x  5 

Yf,  Yg. 
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3.  NUMERICAL  ANALYSIS  OF  THE  PROBLEM 


The  numerical  analysis  of  this  problem  consists  primarily 

of  solving  the  differential  equation  system  (46)  -  (48)  by 

numerical  means  and  secondarily  satisfying  the  edge  conditions. 

We  shall  deal  below  with  these  two  aspects  of  the  problem  and 

then  describe  other  aspects  of  the  computer  program. 

The  treatment  of  the  differential  equation  system  is 

based  on  the  well-known  Runge-Kutta  method  for  numerical 

integration.  The  procedure  is  not  quite  routine  because 

of  the  need  to  evaluate  nc  and  insert  it  in  the  right  side 

of  (47)  and  (48)  at  each  step,  n  being  determined  by  (50)  - 

c 

(52).  In  particular,  if 

l+qay-  ( *y  /  (2x)  }<0  , 

we  must  determine  B  by  solving  (51)  and  then  evaluating  t 

c 

from  (50).  The  solution  of  (51)  was  accomplished  by  using 
the  Newton  -  Raphson  iteration  defined  by 


As  usual  with  this  procedure  convergence  is  uncertain.  However, 
it  became  clear  after  a  little  study  that  convergence  occurred 
when  the  initial  B  was  chosen  well  above  the  value  of  g  for  which 
cos  s  ■  y/x. 

This  condition  was  satisfied  by  choosing  as  starting  values 
B  **  .4  if  c  ^0 
B  -  .4+(-2e)1/2  if  e  <0 . 
e  *  (y  /  x)  -  1 . 
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In  order  to  avoid  difficulties  for  small  values  of  B 


in  (50),  explicit  polynomial  approximations  were  used  for  the 

trigonometric  function,  i.  e., 

2  4 

cos  fi  ^  1-C, 8  +C  8 
1  2 

sin  8  jsSfl-S  82+S  84} 

1  2 

C.  ,  C  ,  S,,  S  constants 
*•  2  1  2 

With  these  approximations  we  may  write  (66)  and  (50)  as 


6k+l  CT  Bk+ 


S  qay-8{  e+82(S1-S2B2)] 


e+82  (C, -C  B2) 
1  2 


z  8k  -  Hk  (67) 


7  2 
(e+BZ(S  -SjB  )) 


U-b2(si-s2b2)  } 


2  4 

{ 1-C  B  +C^  B  ) 


(68) 


These  formulas  are  sufficiently  accurate  for  practical  purposes 
and  are  well-behaved  for  0-$B$it/2. 

Because  the  edge  conditions  (54)  -  (56)  are  not  all  applied 
at  the  same  point,  it  is  necessary  to  adopt  a  trial-and-error 
scheme  in  order  to  satisfy  them.  This  is  begun  by  assuming  a 
value  for  the  deformed  vent  radius,  y(l).  Next  the  values  of 
ns(l)  and  $(1)  are  found  from  (54)  and  (55),  nc  (1)  is  evaluated 
as  in  the  preceding  paragraph  and  we  can  then  begin  the  Runge- 
Kutta  procedure  for  integrating  the  differential  equation  system. 
With  this  procedure  the  solution  is  built  out  in  the  x-direction 
until  the  values  at  the  skirt,  y(x  ),  4>(x  )  and  n  (x  )  are  obtained. 

u  o  So 
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These  values  should  satisfy  (56).  Usually,  of  course,  they  will 
not,  and  we  must  change  the  assumed  value  of  y(l)  and  start  aeain. 
In  the  program  this  process  is  made  automatic  by  including  a 
subroutine  that  applies  the  Rule  of  False  Position  when  a  sign 
change  is  found  in  the  left  side  of  (56). 

In  the  program  a  general  pressure  distribution  is  assumed, 
given  by 


The  constants  A  ,  A. ,  A,  and  A  are  read  in  and  permit  a  fairly 

o  i.  z  3 

wide  choice  of  pressure  distributions.  The  scarcity  of  reliable 
information  on  pressure  distributions  makes  this  procedure 
necessary.  It  is  worth  noticing  that  is  assumed  to  be  a 
function  of  <p  (the  deformed  cord  angle)  and  so  depends  to  a 
certain  extent  on  the  deformed  shape. 

After  finding  the  complete  solution,  satisfying  all  the 
edge  conditions,  the  program  carries  out  the  numerical  integration 
of  (58),  using  Simpson’s  Rule,  and  calculates  the  cylindrical 
coordinates  of  the  cord  and  gore  centerlines  by  means  of  (59)  - 
(62). 

In  order  to  make  the  program  useful  for  design  purposes, 
provision  is  made  for  the  computation  of  certain  information 
beyond  merely  the  shape  and  stresses.  This  extra  information 
includes  an  analysis  of  the  parachute  weight,  i.  e.,the 
distribution  of  the  total  weight  among  the  different  structural 
elements  (fabric,  cords,  fixtures  etc).  Also  the  maximum 
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cord  and  fabric  stresses,  as  well  as  the  vent  and  load  line 
stresses,  are  calculated  and  compared  with  the  respective 
breaking  stresses  to  find  safety  factors  for  each  structural 
element.  Warning  messages  are  printed  if  the  safety  factors 
fall  below  unity. 

The  program  is  arranged  so  that  in  one  pass  the  user  may 
choose  any  two  of  the  input  quantities  and  vary  these 
independently  by  chosen  numbers  of  chosen  increments,  so  that 
the  effect  on  the  solutions  of  changes  in  these  quantities 
may  be  studied.  For  instance,  the  drop  velocity,  V,  and  fabric 
modulus,  E  f ,  may  be  varied,  or  the  effect  of  changing  two  constants 
in  the  hypothetical  pressure  distribution  (69)  may  be  found. 
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4.  EXAMPLES  OF  RESULTS 

This  program  was  used  to  study  the  effects  of  different 
variables  on  the  shape,  drag  and  stresses.  For  this  purpose  a 
basic  set  of  parameter  values  was  chosen,  typical  of  the  C-9  canopy. 
Vent  radius  ■  1.4  ft 
Ratio  of  skirt  to  vent  radii  *  10 
No.  of  gores  ■  28 
Suspension  Line  Length  »  28  ft 

3 

Elastic  Modulus  of  Cords  *■  2  x  10  lbs 

3 

Elastic  Modulus  of  Fabric  •  2  x  10  lbs/ft 
Elastic  Modulus  of  Load  Lines  *  2  x  103  lbs 

The  mass  density  of  air  was  taken  as 
0  -  2  x  10"3  lb-sec/ft4. 
and  the  drop  velocity  as 
V  ■  20  ft/sec 

The  basic  pressure  distribution  was  taken  as  a  constant  function 
C  ■  1.5,  defined  by 

r 

Aq  “  1.5,  Aj  -  Aj  -  A3  -  0 

For  this  basic  configuration  Che  following  results  were  found: 

C  ■  Drag  coefficient  *  .626 
D 

o 

D  •  max.  diameter  »  .671  D 
P  0 

o 

$  ■  cord  angle  at  skirt  ■  109 

or  ° 

9  -(*/2)  -  19° 

o 
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where  D  ■  2R  is  the  flat  circular  diameter  and  ®  is  the  angle 
o  o 

between  the  suspension  lines  and  the  canopy  axis.  The  deformed  shape 
is  shown  in  Figure  13,  and  the  cord  and  fabric  stress  coefficients 
are  depicted  in  Figure  14. 

The  theory  was  further  tested  by  investigating  the  influence 
of  various  parameters  on  the  shape,  drag  and  stresses.  A  large 
number  of  results  was  obtained,  the  most  important  of  which  are 
displayed  in  Tables  1  and  2.  We  may  summarize  the  main  features  as 
f o  1  lows : 

(i)  Table  1  shows  that  D0 /D  ,  and  hence  the  shape,  is  not 

r  o 

greatly  affected  by  any  of  the  parameter  changes.  The  suspension 
line  length  ratio.  A,  has  the  greatest  effect  on  Dp/DQ,  but  its 
effect  is  not  especially  large. 

(ii)  According  to  Table  1,  A  is  the  only  canopy  parameter 
that  affects  the  drag  coefficient  significantly.  The  effects  of 
G  and  V  are  small,  and  the  elastic  constants  have  negligible  influence. 

(iii)  The  angle,  0,  of  the  load  lines,  which  is  the  same  as  the 
angle  of  the  i ords  at  the  skirt,  is  seen  in  Table  1  to  be  significantly 
affected  only  by  changes  in  A- 

(iv)  We  conclude  from  Table  1  that  G  is  the  only  variable  that 

greatly  influences  C  although  A  has  some  effect.  In  contrast 

c  max 

Cf  max  is  perceptibly  affected  by  changes  in  all  the  canopy  parameters. 

(v)  The  cord  stress  is  least  at  the  vent  and  increases  to  a 

maximum  at  the  maximum  diameter,  which  is  usually  at  x  *  9,  i.  e. 

near  the  skirt.  Hence  the  familiar  formula 

Cc  "  CD  /G  cos  0, 
o 
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E 

r 


TABLE  2 


In  the  case  labelled  "Linearly  varying  C  " ,  C  varies  linearly  In  $ 
between  Cp ( 1  }  -  1.5  (at  vert)  and  Cp(xQ)Pat  tSe  skirt. 


which  is  exactly  true  only  for  C  (x  )  ,  it  also  approximate ly  true 

c  o 

for  the  maximum  value  of  C  . 

c 

(vi)  The  fabric  stress  is  greatest  at  or  near  the  vent  and 
decreases  as  we  go  toward  the  skirt.  The  exact  location  of  the 


maximum  fabric  stress  dependson  E ^ 
diet  ribution  • 


as  well  as  the  pressure 


(vli)  From  Table  2  we  see  that  for  constant  pressure  distributions. 


changing  the  pressure  scarcely  affects  the  shape  and  causes  merely 


proportional  changes  in  the  stresses.  If  the  pressure  Increases 
linearly  in  <P  from  the  vent  to  the  skirt,  all  quantities  are  affected. 


Ue  shall  comment  on  these  results  in  the  following  Section. 
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5.  DISCUSSION 


The  theory  developed  in  this  paper  is  apparently  the  first 
in  which  the  shape  and  stresses  are  calculated  simultaneously. 

To  Judge  the  theory  we  must  of  course  compare  it  with  experimental 
results.  In  making  this  comparison  it  is  first  necessary  to 
understand  clearly  what  the  theory  does. 

If  we  are  given  a  flat,  circular  canopy  with  known  geometrical 

and  physical  properties  (respectively  embodied  in  R.  ,  R  ,  L  ,  0 

i  o  0 

and  E^,  Ec,  and  E^) ,  dropping  vertically  at  known  velocity  V  through 
air  with  known  mass  density,  p,  the  analysis  furnishes  us  with  the 
relations  between  the  oressure  distribution  on  one  hand  and  the 
shape*  stresses  and  drag  on  the  other.  If  the  pressure  distribution 
is  known,  we  can  solve  these  relations  (by  means  of  the  computer 
program)  to  find  the  quantities  of  engineering  interest.  These 
quantities  will  depend  in  general  upon  the  pressure  distribution. 

An  ideal  experimental  check  on  the  theory  would  require  that 
simultaneous  measurements  of  pressure  distribution,  stresses,  shaoe 
and  drag  be  made  on  a  canopy  that  is  dropping  vertically.  Conditions 
for  such  a  test  may  be  difficult  or  impossible  to  realize  in  practice, 
and  it  appears  that  no  experiments  yet  made  will  permit  such  a 
complete  check  on  the  theory. 

However,  many  experiments  have  been  made  that,  while  incomplete 
in  some  respect,  give  us  the  information  for  a  partial  check  on  the 

3 

theory.  For  example,  Berndt  and  Deweese  made  measurements  on  a 
C-9  canopy  in  towed  flight.  Their  results  compare  with  the  present 
calculation  for  the  basic  configuration  as  follows: 


Meas  ured 


Calculated 


C 


D 


.65 


o 


C  -  .626 
D  _ 


D  /D  -  .648 
P  o 


D  /D  -  .671 
P  ° 


-  17* 


0  -  19° 


Their  measured  shape  is  plotted  in  Figure  13.  The  general  agreement 

Is  good  although  the  theoretical  shape  is  slightly  wider  and  shallower 

than  the  experimental  shape.  We  may  conclude  from  this  that  the 

present  theory  is  not  grossly  in  error.  However,  It  Is  necessary 

to  remark  that.  If  we  had  assumed  a  different  pressure  distribution 

In  the  calculations,  we  could  have  arrived  at  quite  different  values 

for  C  ,  as  is  evident  from  Table  2. 

D 

°  4 

The  Parachute  Handbook  gives  the  general  estimate 

D  /D  «  .7, 

P  ° 

which  compares  reasonably  well  with  the  present  calculations, 

D_/D  ■  .67  for  the  cords 
P  o 

*  .73  for  the  gore  centerlines. 

Moreover,  the  observed  generality  of  this  estimate  tends  to  confirm 
the  theory,  which  predicts  that  the  shape  Is  Insensitive  to  changes 
in  most  of  the  parameters. 

4 

The  Parachute  Handbook  shows  many  experimental  results  In  which 
C  Is  seen  to  decrease  with  an  Increase  in  V.  This  does  not  agree, 

n 

at°firsr  sight,  with  our  results.  Table  1,  which  show  that  Cn 

uo 

increases  very  slowly  with  an  increase  in  V.  The  discrepancy  is 
probably  due  to  one  or  both  of  two  causes: 
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(a)  The  canopy  may  not  be  descending  vertically  in  the  tests. 

The  theory  assumes  vertical  descent  and  may  give  inaccurate  results 
if  the  canopy  oscillates  or  glides  laterally. 

(b)  The  pressure  distribution  may  change  as  the  velocity 

changes.  In  particular,  increasing  velocity  will  cause  increase 

of  the  tensile  stresses  and  strains  in  the  fabric,  and  this  will 

lead  to  greater  porosity  of  the  fabric.  The  increasing  norosity 

will  permit  increasing  air  flow  from  the  high-pressure  to  the  low- 

pressure  side  of  the  fabric.  Thus  the  pressure  difference  will  not 

increase  as  fast  as  the  impact  pressure,  or,  in  other  words,  the 

pressure-difference  coefficient,  Cp ,  will  decrease.  This  will 

cause  to  decrease, 

o 

In  general  we  may  say  that  the  present  theory  represents  a 
first,  admittedly  somewhat  crude,  attempt  to  analyze  parachutes 
without  assuming  the  deformed  shape.  Clearly,  several  of  the 
assumptions  described  in  the  Introduction  are  not  very  accurate, 
notably  (ii),  and  (v)  ,  while  (iv)  is  probably  inaccurate  near  the 
maximum  diameter.  Nevertheless  the  results  are  good  enough  to 
suggest  that  the  general  procedure  is  a  workable  one  which  can  be  a 
starting  point  for  future  efforts. 

It  is  perhaps  also  useful  to  point  out  that  this  theorv  represent 
a  step  in  the  direction  of  calculating  the  opening  behavior  of  para- 
chutes.  For  a  major  part  of  the  difficulty  in  the  opening  problem 
consists  of  estimating  the  shape  and  forces  of  the  inflated  part  of 
the  canopy.  The  no t -ye t -in f lat e d  part  of  the  canopy  is  relatively 
easily  dealt  with.  The  estimate  of  the  shape  and  forces  of  the 


Inflated  part  of  the  canopy  is  much  like  that  of  a  (smaller) 
completely  open  canopy.  We  may  anticipate,  therefore,  that  the 
present  analysis,  or  one  like  it,  will  be  an  integral  part  of  an 
analysis  of  opening  behavior. 
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NOMENCLATURE 

Constants  In  canopy  pressure  distribution 
Constants  in  formulas  for  cosine  function 

Drag  coefficient  based  on  flat  circular  area 

Coefficient  of  net  pressure 

Stress  coefficients  in  cords,  fabric 

Diameters  in  opened  and  flat  circular  states 

Elastic  moduli  of  cords,  fabric,  load  lines 

1+<N  /E) 
c  c 

Number  of  gores 

Dimensionless  load  line  length 

Length  of  load  lines  in  undeformed  and 
deformed  states 

Tension  force  in  cords,  fabric 

Force  resultants  on  cords  in  circumferential 
and  normal  directions 

Dimensionless  force  in  cord,  fabric 

Dimensionless  circumferential  force  resultant 
on  cords 

Net  outward  pressure  on  gore  fabric 
Dimensionless  pressure 

Radius  of  particle  *n  flat  circular  state 
Vent  and  skirt  radii  in  flat  circular  state 
Radius  of  cord  in  opened  state 
Vent  and  skirt  radii  in  opened  state 
Circumferential  radius  of  curvature  of  gore 
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NOMENCLATURE  (Cont’d) 


s 


S 


1 


S 

2 


Arc  length  of  cord 

Constants  In  formulas  for  sine  function 


Dimensionless  elastic  constant* 

Vertical  velocity 
Drag  force 
Dimensionless  drag 

Dimensionless  radius  In  flat,  circular  state 

Dimensionless  skirt  radius  In  flat  circular 
state 


y 


Dimensionless  radius  in  opened  state 


7., 


Axial  cylindrical  coordinate  of  cords  and 
gore  centerline 


z  . 


x 


R 


Dimensionless  axial  coordinate  of  cords  and 
gore  centerline 


Included  angle  of  gore 

Edge  angle  of  gore  bulge 

Strains  in  cord,  fabric,  and  load  lines 

Depth  In  bulge 

Dimensionless  depth  of  bulge 

Angle  between  load  line  and  axis 

Load  line  length  ratio  -  T,  /D  -  JR-/(?R  ) 

o  o  1  o 

Mass  density  of  air 

Contact  length  between  adjacent  gores 

Dimensionless  contact  length  between  adjacent 
gores 

Angle  between  cord  direction  and  horizontal 
$  at  the  skirt 
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